Electric-magnetic duality and topological order on the lattice 
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Dualities are deep conceptual tools in many areas of physics and mathematics; in particular, 
electric-magnetic (EM) duality has a long tradition ranging from Dirac's magnetic poles to S-duality 
in string and field theory. Here we investigate the duality structure in quantum lattice systems 
with topological order, a collective order also appearing in fractional quantum Hall systems. We 
define EM duality for all of Kitaev's quantum double models based on discrete gauge theories with 
Abelian and non-Abelian groups, and identify its natural habitat as a new class of topological models 
based on Hopf algebras. We interpret these as extended string-net models, whereupon Levin and 
Wen's string-nets, which describe all topological orders on the lattice with parity and time-reversal 
invariance, arise as magnetic and electric projections of the extended models. We conjecture that 
all string-net models can be extended in an analogous way, using a more general algebraic structure, 
such that EM duality continues to hold. Physical applications include topology probes in the form 
of pairs of dual tensor networks. 



Duality ranks among the deepest ideas in physics. Du- 
alities are powerful probes into the internal structure of 
mathematical and physical theories, including quantum 
many-body systems and field theories [1]. Frequently a 
duality relates weakly and strongly coupled regimes of 
physical systems, providing precious information beyond 
the reach of perturbative methods. Many physical in- 
stances of duality involve the exchange of electric and 
magnetic degrees of freedom. Such a symmetry for the 
Maxwell equations in vacuo led Dirac [2] to the intro- 
duction of pointlike sources of magnetic field to extend 
this electric-magnetic (EM) duality to matter, providing 
a unique argument for the quantisation of electric charge. 
Dirac's insight has had profound influence in field theory 
and beyond (see, e. g., [3]). 

Topological order [4] appears in a fascinating class of 
condensed matter phases, notably the fractional quan- 
tum Hall effect. Topological phases are sensitive only 
to global properties of the underlying space: their effec- 
tive quantum field theories have only non-local degrees of 
freedom [5] . This non-local order has deep connections to 
high-energy physics and a rich mathematical structure. 
In 2D, it is associated with localised objects with frac- 
tional statistics, that is to say, anyons [6]. In topological 
lattice models, these anyons appear as excitations. 

Kitaev's proposal [7] to use topologically ordered quan- 
tum lattice systems to store and process quantum in- 
formation linked topological order and quantum com- 
putation. His pioneering quantum double D(G) models 
are built in analogy to gauge theories with finite gauge 
groups. In this context, the topological degrees of free- 
dom are Wilson loops; violations of gauge invariance and 
magnetic fluxes are described analogously, hinting at EM 
duality. Particularly simple are the models based on 



Abelian groups, including the toric code, the first labora- 
tory for many ideas in topological order. Although there 
are deeper objects at play, D(G) models can be conve- 
niently studied using group theory; there is an explicit 
local understanding of anyonic excitations and their in- 
terplay in terms of group representations. These models 
provide us with the first emergence of duality in topolog- 
ically ordered systems, namely a self-duality of the toric 
code and Abelian D(G) models, involving the direct and 
dual lattices; an EM duality, in gauge theory language. 
Its importance has been stressed by Fendley [9] . 

Levin and Wen defined their landmark string-net (SN) 
models [8] from the physical intuition that topological 
phases are renormalisation group fixed points. SNs are 
built on ideas from category theory to describe all topo- 
logical orders on the lattice with parity and time-reversal 
invariance. Technically, they use data from unitary ten- 
sor categories; in this spirit, dualities in SN models, in 
the sense of equivalences of unitary modular tensor cate- 
gories, can be identified with invcrtiblc domain walls [ Hi]. 
While SN models are broader in scope than Kitaev's 
quantum double models, their excitations are, however, 
not understood as thoroughly as those of the D(G) mod- 
els, especially locally. 

So far the extension of the EM duality of the toric 
code to non-Abelian models has remained an open prob- 
lem, and thus a key source of insight into the structure 
of topological order has been only available in the very 
restricted Abelian setting. The natural place to look for 
a well defined EM duality are Kitaev's quantum double 
models, since the excitations there are clearly understood 
and can be assigned electric and magnetic quantum num- 
bers, while this characterisation is absent in SN models. 

Here we uncover EM duality for all of Kitaev's D(G)'s, 
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FIG. 1. A duality landscape. We represent: Quantum dou- 
bles T)(A) based on Abelian groups, including the toric code, 
all of which are self-dual; general group QD models D(G), 
related by duality to quantum double models D(G*) based 
on algebras of functions on groups, whose intersection with 
D(G) is D(j4); these are all instances of quantum double mod- 
els D(_ff) based on Hopf algebras, a larger class closed under 
duality, with cases of self-duality beyond groups [11]. We con- 
jecture that EM duality can be defined for extended string-net 
models ESN, and that these can be identified as quantum dou- 
ble models D(W) based on weak Hopf C*-algebras. Arrows 
denote model identifications upon lattice dualisation. 



both Abelian and non- Abelian, and beyond. We do this 
by identifying its natural setting as the quantum double 
models based on Hopf algebras, a class of models an- 
ticipated in [7] and defined and studied in [11]. Then 
we use Fourier techniques [12] to write these as extended 
string- net models (ESNs), a variant of Levin and Wen's 
SNs with a canonically enlarged Hilbert space, such that 
anyons can be understood in representation theoretic lan- 
guage [27]. ESNs can be defined independently, in anal- 
ogy to Levin and Wen's SNs, by data from unitary tensor 
categories, equipped in this case with a so-called fibre 
functor [13]. By EM duality, each ESN maps to two 
SNs, its electric and magnetic projections (as depicted 
in figure 6). We conjecture that all SNs are obtained in 
such a way from parent ESNs with a generalised alge- 
braic structure, in which EM duality continues to hold. 
These results lead us to propose the duality landscape of 
figure 1. 



KITAEV'S MODELS BASED ON GROUPS 

Let us begin by reviewing the well-known self-duality 
in the toric code [7] , the quantum double model based on 
the Z2 group. This is a model of qubits along the bonds of 
a 2D lattice A (the dimension of the local Hilbert spaces 



being the size of the group), with a Hamiltonian of the 
form 



H 



(i) 



where v runs over vertices and p runs over plaquettes of 
A. Mutually commuting vertex and plaquette operators 
are projectors involving Pauli matrices: 



A v = 



^ , £>?j — ^ 



, (2) 



with support on the bonds around the corresponding ver- 
tex or plaquette. Self-duality means that, upon a global 
unitary 



u A = 00 u i . & = 



ieA 



(<J x + <J z )i 
V2 



(3) 



built from Hadamard maps on each edge i, the toric code 
on A gets mapped to a toric code on the dual lattice A* : 



U A A v Ul = B v , U A B p U{ = A p 



(4) 



where now B v corresponds to dual plaquette v, and A p 
to dual vertex p in A*. Ground states of the frustration- 
free Hamiltonian (1) minimise each term in the sum, 
that is, they satisfy all vertex and plaquette constraints 
A v \tp) = \ip), B p \tp) = \tp). Breakdown of any one such 
constraint means the state is in an excited level, and can 
be interpreted as the presence of a localised excitation, or 
particle, sitting at the vertex or plaquette that is broken; 
from the gauge theory interpretation, these are called 
electric and magnetic excitations. The electric-magnetic 
nature of duality (4) comes then from the interchange of 
vertices and plaquettes in going from the original to the 
dual lattice. 

Hidden in the action of the global Hadamard map Ua 
arc, first, the mapping from A to the dual lattice and 
the corresponding reinterpretation of vertex and plaque- 
tte operators, and second, a mapping of elements of the 
group Z2 to the functions from this group to the scalars. 
If we obtain again a toric code in the dual lattice it is 
because, for an Abelian group, the space of functions has 
the same structure as the group algebra, which is the 
space of linear combinations of group elements. One can 
exploit this algebraic fact to extend the self-duality of 
the toric code to all D(G) models for Abelian groups in a 
straightforward way: this is the D(A) region of figure 1. 

For general D(G) models, we start with a lattice with 
oriented bonds, whose Hilbert spaces have an orthonor- 
mal basis {\g)} labelled by elements of a finite group 
G. The Hamiltonian of the D(G) model is still of the 
form (1), but now the mutually commuting vertex and 
plaquette operators are the projectors defined in figure 2. 
The obvious parallelism of plaquettes and vertices in the 
Abelian case is lost if G is non- Abelian. Upon switching 
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FIG. 2. Projectors in the D(G) model. Vertex (A v ) 
and plaquette (B p ) projectors in the Hamiltonian of Kitaev's 
D(G) model are defined by their action on a basis of the 
corresponding multi-edge Hilbert space. Oriented edges are 
labelled by group elements denoting kets in the orthonor- 
mal computational basis {g} a eC, change of orientation is im- 
plemented by applying the inversion map g i-> g 1 - The 
vertex operator projects onto gauge- invariant configurations 
(=symmetrises over simultaneous multiplication with group 
elements). The plaquette operator involves a Kronecker delta 
on the unit element e £ G and projects onto trivial magnetic 
flux across p. 



to the dual lattice, we cannot reinterpret the operators 
A v and B p as plaquette and vertex operators in a quan- 
tum double model based on a group, since the functions 
on a non-Abclian group cease to have the same structure 
as the group algebra itself. Hence, the extension of the 
EM duality to non-Abelian groups is not possible within 
the class of quantum doubles based on groups. Yet the 
construction of the Abelian duality can be generalised 
to the non-Abelian case, resulting in dual models which 
are no longer quantum double models based on groups, 
but rather quantum double models based on algebras of 
functions; this is shown as the arrow between the regions 
D(G) and D(G*) in figure 1. To make sense of this, one 
has to widen the construction of quantum double lattice 
models beyond the group case. 

As we will show, the natural habitat for the EM duality 
of D(G) models is the class of quantum double models 
based on Hopf algebras [11]; this is the region D(£f) in 
figure 1. This is the smallest setting containing all D(G) 
models and closed under tensor products and EM duality; 
the latter, moreover, takes a remarkably simple form in 
the language of Hopf algebras. 



THE D(H) LATTICE MODELS 

When studying quantum many-body systems, more 
general notions of symmetry emerge than those furnished 
by group theory. In particular, linear transformations 
acting on tensor products of vector spaces lead naturally 
to Hopf algebras. For a detailed account of these in our 
context wc refer to [11]. In the following we give an 



intuitive grasp of their structure, which is necessary to 
understand the D(H) lattice models. 

We consider Hopf algebras H as spaces of transforma- 
tions on a many-body Hilbert space. First of all, the 
linear nature of the target is naturally extended to its 
transformations, so Hopf algebras are vector spaces. We 
must be able to compose transformations and to include 
the identity transformation, so H has a multiplication of 
vectors, and a unit. Most importantly, we must have a 
rule to distribute the action of an element of H into a 
tensor product of target spaces; this is the so-called co- 
multiplication. Additionally H has a trivial representa- 
tion e, called counit, making precise the notion of spaces 
invariant under the action of H. Like for groups, the rep- 
resentation theory of Hopf algebras includes a notion of 
conjugate representation, implemented via an antipode 
mapping, which for groups is just the inversion g i— > g^ 1 . 

In order to construct Hilbert spaces, we use Hopf C*- 
algebras, which in addition have an inner product. We 
will only consider finite-dimensional Hopf G*-algebras, 
which come equipped with a canonical, highly symmet- 
ric element, the Haar integral h € H. invariant under 
multiplication in the sense a ■ h = e(a) h for all elements 
a of H. This canonical element is crucial for the con- 
struction of the lattice model and its ground states. 

The root of the EM duality to be unveiled in the fol- 
lowing is an algebraic fact: The class of Hopf C* -algebras 
is closed under dualisation. That is, given a Hopf G*- 
algebra H, its dual space H* of functions from H to 
the scalars is again a Hopf G*-algebra, whose structure 
is, moreover, determined by the structure of H, as illus- 
trated by figure 3. This closure property is shared by the 
class of Abelian group algebras (which are all self-dual) , 
but not for the whole class of group algebras. The land- 
scape of figure 1 reflects these statements in the world of 
lattice models. 

Figure 4 defines a set of operators on a lattice whose 
bonds carry elements of H , a finite-dimensional Hopf G*- 
algebra. These operators depend on elements of the al- 
gebra and its dual, respectively. When these are taken 
as the Haar integrals h of H and </> of H*, the resulting 
operators 
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A ViP (h) , B p = B VtV (4>) 



(5) 



define via (1) the Hamiltonian of a topological lattice 
model: the quantum double model based on H, or D(iJ) 
model for short [11]. 

Kitaev's models are a particular subclass of D(i/) mod- 
els, recovered when H is the group algebra of a finite 
group G. Both this algebra CG and its dual G* have par- 
ticularly simple Hopf G*-algebra structures, summarised 
in figure 3. In particular, their Haar integrals read 
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FIG. 3. Structure and dualities in Hopf algebras. The 

first row represents the structure maps of a Hopf algebra H. 
Solid lines represent elements of H, and dotted lines repre- 
sent scalars in C. The second and third rows illustrate the 
algebraic duality within the class of Hopf algebras, taking as 
examples a group algebra H — CG and its dual H* , whose 
structure is defined implicitly in the figure. For each dia- 
gram, the upper expression is equal to the lower expression; 
the structure maps of H and H* are represented by solid and 
dashed lines, respectively, according to the conventions of the 
first row. Hence, the multiplication in H* is determined by 
the comultiplication in H, the unit in H* is determined by 
the counit in H* , etc.; this duality holds for general finite- 
dimensional Hopf algebras. We use the form of the structure 
maps in the basis {g; g £ G} of the group algebra, and the ba- 
sis {8f, I £ G} of its dual, where St : g >-¥ 5i, g are Kronecker 
delta functions. As for the extra structure needed to define 
our lattice models, the C*-algebra structures of H and H* 
are defined by declaring {g} to be an orthonormal basis, and 
{Si} to be orthonormal up to a factor, (Si, S m ) = |G| _1 <^, m . 
The Haar integrals h and (j> in the group algebra and its dual 
are defined in equation (6). 



Taking into account this equation, operators in figure 2, 
defining the Hamiltonian of the D(G) model, can be 
checked to coincide with the general expression (5). 
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FIG. 4. The core of D(ff) models. Topological models 
based on Hopf algebras are defined via these vertex and pla- 
quette representations of u £ H and / £ H* [28]. Oriented 
edges carry, in general, elements of a finite-dimensional Hopf 
C*-algebra H; for Kitaev's models, take H to be the group 
algebra CG. Orientation can be reversed by using the an- 
tipode map S referred to in the text. The sums denote the 
appropriate comultiplication or splitting of algebra elements 
into several tensor factors. Vertex operators A VtP now de- 
pend also on a plaquette p marking the beginning of a loop 
around the vertex v, and correspondingly for the plaquettes. 
These operators probe the algebraic structure of the mod- 
els: Given v and p, A VlP (u) and B p . v ( f) interact nontrivially, 
building representations of Drinfel'd's quantum double algebra 
D(_ff) [14]. The irreducible representations of D(_ff) classify 
the superselection sectors, and thus the anyonic structure of 
excitations, of the quantum double lattice models [7] . 



ELECTRIC-MAGNETIC DUALITY 

We now define EM duality for general D(H) models. 
Consider the following unitary map Li from H to its dual: 



la , f a (b) 



\H\4>(a-b) , 



(7) 



where \H\ is the dimension of H, and the function f a 
is constructed via the dual Haar integral (j> £ H* (for 
instance, for groups, f g = y/\G\ S g -i). We associate this 
mapping with the transformation of A into its dual lattice 
A* as shown in figure 5. 

The vertex and plaquette representations in figure 4, 
associated with the D(if) model on A, are mapped by 
= ® ieA ^i into precisely the representations associ- 
ated with the B(H*) model on the dual lattice A* [29]: 

U A B p . v {f)u{ = Ap iV (f) , U A A VtP (u)Ul = B v . p {u), 

(8) 

where p is a plaquette of A and a vertex of A*; this ex- 
pression generalises equation (4) and identifies Ua as the 
EM duality mapping. 

Thus, U\ identifies the D(_ff) model on A with the 
D(H*) model on A*, as it transforms the Hilbert spaces 
and Hamiltonian of the former into those of the latter. 
We write this symbolically as [D(H)] A = [D(i3"*)] A «. 

But EM duality is deeper than a mere identification of 
energy levels. The Hamiltonian involves just the opera- 
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FIG. 5. Going over to the dual lattice. The mapping 
U : H — > //* is here interpreted as a lattice dualisation. The 
solid bond is in A, while the dashed bond is its dual in A*. 
The orientation of the dual edge is chosen by convention so 
that it points to the left of the direct edge. The degree of 
freedom a £ H in the Hilbert space H of the bond in A 
becomes a degree of freedom /„ of the Hilbert space H* of 
the dual bond. When applied simultaneously on all edges of 
A, this mapping implements EM duality. 

tors A v (h) and B p (<p), but, according to (8), Ua induces 
a transformation of all operators of figure 4 that respects 
their algebra structure, and relates the anyonic excita- 
tions of both models [13]. 

FROM QUANTUM DOUBLES TO 
STRING-NETS 

We now make the connection with Levin and Wen's 
string-net models [8]. In [12] the D(G) model was identi- 
fied, via a Fourier mapping, with an extended string-net 
model on the same lattice, with bond degrees of free- 
dom given by irreducible representations (irreps) of G 
and auxiliary matrix indices depending on each irrcp (sec 
also [15]). String-net vertex conditions amount to knit- 
ting the matrix indices together by means of 3j symbols, 
effectively mapping the local Hilbert spaces isometrically 
to those of a Levin- Wen SN model, whose degrees of free- 
dom and underlying category theoretical structure are 
given solely by irreps of G and their fusion properties. 
The ground levels of SN and their extended versions are 
identical. 

Now the representation theory of finite Hopf C*- 
algebras is essentially identical to that of finite groups 
(the only real difference being that the fusion of rep- 
resentations need not be commutative). Therefore the 
construction of [12] can be generalised immediately to 
any D(7J) model. That is, if H has matrix irreps of 
dimension we define a basis in H as 

lt,ab] = V\HjtiJ2 D a b ( hil) ) h{2) , (9) 

CO 

where a, b = 1, . . . , \fj,\ are matrix indices for irrep D^. 
This is an orthonormal basis for each bond, and the 
analysis of [12] carries through intact to show that the 
[D(J?)]a model can be written as an extended string- 
net model [ESNh-]a- The latter reduces to a string- net 



model [SN^ ag ]A, whose degrees of freedom are the irreps 
of H, and whose ground level is identical to that of the 
extended model. 

The same construction, as applied to the [D(H*)]\* 
model using the Fourier basis for H* , yields an equivalent 
extended string-net [ESN#.]a* on the dual lattice. This, 
in turn, reduces to another string-net model [SN^]a* on 
the dual lattice, whose bond degrees of freedom are irreps 
of H* . We regard the two string-net models [SN^ ag ]A 
and [SN^]a* as the magnetic and electric projections of 
the original quantum double model (see figure 6). By 
construction, they satisfy [SN^a* = [SN™ ag ] A *. 

For instance, in [12] the magnetic SN projection of the 
D(G) model was analysed. Its electric SN projection, on 
the other hand, has the same local degrees of freedom 
as the D(G) model, since group elements are the (one- 
dimensional) irreps of the dual of the group algebra [30] . 

Thus, EM duality can be given a meaning in the sub- 
class of string-nets obtained by reduction of ESNs. This 
leads us to the following 

Conjecture: All SN models can be extended into quan- 
tum double models based on a sufficiently general al- 
gebraic structure, the class of weak Hopf G*-algebras 
[16] [31]. EM duality can be defined for all extended SN 
models so that the pattern of figure 6 holds. 

PROBING TOPOLOGY WITH TENSOR 
NETWORKS 

Some physical consequences of EM duality are imme- 
diate. First of all, EM duality allows us to probe the 
topology of the surface underlying the lattice models by 
using only locally defined states. 

As shown in [11], each [D(H)}\ model has one canoni- 
cal tensor network ground state \i/j(H, A)) constructed 
from identical tensors at each site, which are defined 
solely by the structure of H. 

From the corresponding canonical state \ip(H*, A*)) of 
the dual [D(.H"*)]a* model on the dual lattice we can ob- 
tain another ground state of the original [D(H)]\ model 
by using the duality map (7): 

A)) = U A \ip(H*, A*)) . (10) 

The relation between \4>(H, A)) and \ip(H, A)) de- 
pends on the topology of the surface underlying A. On 
the sphere, for instance, the ground level is nondegener- 
ate, so both ground states are the same. On a topolog- 
ically nontrivial surface such as the torus, these ground 
states can be shown to always be linearly independent. 

For instance, in the toric code the tensor network con- 
struction of the canonical state \ip(H, A)) coincides with 
the projected entangled-pair state (PEPS) ansatz in [18]. 
On the torus this yields the logical 1++) state; the dual 
ground state \ip(H, A)) is the logical |00) state. See [13] 
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FIG. 6. From quantum doubles to string-nets. Quan- 
tum double models and string-net models are interconnected 
by EM duality, the Fourier mapping, and a projection to SN 
degrees of freedom. The direct lattice is drawn solid, the 
dual lattice dashed. The first layer represents the quantum 
double models in both their EM-dual formulations. Via the 
Fourier mapping we write the quantum doubles as extended 
string-net models, depicted in the second layer, which are en- 
tirely equivalent. The degrees of freedom of the ESN model 
on A are labelled by irreducible representations \i of H and 
auxiliary indices a, b; for the dual model on A* the labels are 
irreducible representations Jl of H* , and auxiliary indices a, b. 
These models support excitations both at vertices (electric) 
and at plaquettes (magnetic), which we represent by blobs 
and lozenges. An electric excitation e becomes a magnetic 
excitation fh in the dual, and viceversa. The third layer rep- 
resents the SN models obtained, respectively, by projecting to 
degrees of freedom /x on the direct lattice, or /x* on the dual 
lattice; these we refer to as the magnetic and electric projec- 
tions of the ESN. Excitations in SN models are less symmetric 
than in the ESN: for instance, violation of a vertex constraint 
implies the violation of all surrounding plaquette constraints 
(see [12]). 



for a detailed construction, including the non-Abelian 
cases. 

On the other hand, the correspondence between string- 
net and quantum double models can be used to relate 
these PEPS to the tensor network descriptions of string- 
net ground states put forward in [20] and [21]. These 
come from the \ip(H, A)) construction of the correspond- 



ing [D(.H")]a model, while the \i/j(H, A)) can be seen to 
yield a new TN given simply by the knitting together of 
3j symbols at vertices in the ESN degrees of freedom. 
More details will be given in [13]. The Fourier construc- 
tion is also expected to relate the entanglement renor- 
malization analyses of D(Zf) [11, 22] and SN models [23]. 

DISCUSSION AND OUTLOOK 

We have defined EM duality for non-Abelian D(G) 
models, showing how it arises naturally in the context of 
D(-ff) models. The connection to SN models comes from 
the Fourier construction of extended string-nets and their 
electric and magnetic projections. We conjecture that all 
SN models have extended versions, and that EM duality 
can be defined for all the extended models. EM duality 
offers nonlocal information about the systems, unveiling 
global characteristics of space, from tensor networks de- 
fined locally. Beyond this, it serves as an organising prin- 
ciple for topological models; indeed, we envisage a net of 
dualities for topological phases (not unlike that connect- 
ing theories of strings and branes). This then reflects on 
the field theories underlying topological systems and on 
their mathematical structures, such as category theory. 

We wish to stress the importance of the whole class 
of models based on Hopf algebras. This is indeed the 
natural class of models to study to understand electric- 
magnetic duality for lattice models based on non-Abelian 
groups. We thus shift the focus away from lattice gauge 
theories based on groups, which are a primary concern in 
high energy theory, even when the Hopf algebra language 
is applied [24]. Ours is, moreover, the first proposal for a 
general EM duality governing topologically ordered sys- 
tems. 

The Hopf algebraic language makes the EM duality 
transparent and easy to grasp, instead of obscuring it; in 
addition, this language can be learnt here hands-on from 
the standpoint of quantum many-body systems, which 
can even be built in a laboratory. We also provide, with 
the Fourier construction, a way to understand some as- 
pects of category theory in terms of representations, more 
familiar to many physicists. 

Our next goal is defining the extended string-nets for 
all SN models [13], casting ESNs and their dualities in 
the language of categories. The nature of the anyonic 
models controlling the excitations of the D(i?) models 
has a bearing on the experimental realisation of string- 
net models, which we contend would be most natural in 
their extended incarnations (protocols for anyonic ma- 
nipulation in the spirit of [25] and [26] would then be 
possible). A key question to be explored is the meaning 
of EM upon perturbation of the models considered, that 
is, away from topological fixed points. 

Acknowledgements: We would like to thank 
J. M. Mombclli, J. Slingerland, and G. Vidal for dis- 



7 



cussions. LK is supported by the Gordon and Betty 
Moore Foundation through CPI in Caltech and NSF 
under Grant No. PHY-0803371 and the Basic Research 
Young Scholars Programs of Tsinghua University. 



[1] R. Savit, Rev. Mod. Phys. 52, 453 (1980). 

[2] P. A. M. Dirac, Proc. Roy. Soc. London A133, 60 (1931). 

[3] N. Seiberg, Nucl. Phys. B435, 129 (1995). 

[4] X.-G. Wen, Adv. Phys. 44, 405 (1995). 

[5] E. Witten, Commun. Math. Phys. 117, 353 (1988.) 

[6] F. Wilczek, Phys. Rev. Lett. 48, 1144 (1982); 

Phys. Rev. Lett. 49, 957 (1982). 
[7] A. Yu. Kitaev, Annals Phys. 303, 2 (2003). 
[8] M. A. Levin and X.-G. Wen, Phys. Rev. B71, 045110 

(2005). 

[9] P. Fendley, Annals Phys. 323, 3113 (2008). 
[10] A. Yu. Kitaev and L. Kong, Models for gapped edges and 

boundaries between phases, in preparation. 
[11] O. Buerschaper, J. M. Mombelli, M. Christandl and 

M. Aguado, in preparation. 
[12] O. Buerschaper and M. Aguado, Phys. Rev. B80, 155136 

(2009). 

[13] O. Buerschaper et al., in preparation. 

[14] V. G. Drinfel'd, Quantum groups, in Proceedings of 
the International Congress of Mathematicians, Berkeley, 
California, August 3-11, 1986, ed. A. M. Gleason, AMS, 
1988. 

[15] Z. Kadar, A. Marzuoli, M. Rasetti, preprint 

arXiv: 0907. 3724 [quant -ph] . 
[16] T. Hayashi, preprint arXiv: math/9904073. 
[17] G. Mack, V. Schomerus, Nucl. Phys. B370, 185 (1992). 
[18] F. Verstraete, M. M. Wolf, D. Perez- Garcia, and 

J. I. Cirac, Phys. Rev. Lett. 96, 220601 (2006). 
[19] F. Verstraete and J. I. Cirac, Phys. Rev. A70, 



060302 (2004); F. Verstraete, J. I. Cirac, preprint 

arXiv : cond-mat/0407066. 
[20] O. Buerschaper, M. Aguado and G. Vidal, 

Phys. Rev. B79, 085119, (2009). 
[21] Zh.-Ch. Gu, M. Levin, B. Swingle and X.-G. Wen, 

Phys. Rev. B79, 085118 (2009). 
[22] M. Aguado and G. Vidal, Phys. Rev. Lett. 100, 070404 

(2008). 

[23] R. Koenig, B. W. Reichardt and G. Vidal, 
Phys. Rev. B79, 195123 (2009). 

[24] R. Oeckl, Discrete gauge theory. From lattices to TQFT, 
Imperial College Press, 2005. 

[25] M. Aguado, G. K. Brennen, F. Verstraete and J. I. Cirac, 
Phys. Rev. Lett. 101, 260501 (2008). 

[26] G. K. Brennen, M. Aguado and J. I. Cirac, New 
J. Phys. 11, 053009 (2009). 

[27] The fact that anyon models in string-net models corre- 
spond to representations of algebraic structures was an- 
ticipated in Levin and Wen's original paper [8]. 

[28] Technically, the A v , p (u) axe representations of H, and 
the B PtV (f) are representations of a variant of the dual; 
namely, the dual (H op )* of the opposite algebra of H 
(where the multiplication is flipped). 

[29] Technically, representations A VtP (u) of H and B p , v (f) 
of (H op )* act on bonds of A, that is, spaces of the form 
H® n . They get mapped by Ua to representations B v , p (u) 
of {{H*) op Y and A p , v (f) of H* on bonds of A*, that 
is, spaces of the form (H*)® n . The latter, for coinciding 
v and p, define representations of the D(_ff*) algebra, 
building the D(77*) model on A*. 

[30] This connection of SN and D(G) models had been recog- 
nised long ago by Hector Bombm. 

[31] The excitations in such models would be classified by 
representations of the double D(W / ) of the weak Hopf C- 
algebra W , which is expected to be a weak quasi-Hopf 
algebra in the sense of [17]. 



